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Abstract 

This paper presents some parallel developments in Quiver/Dimer 
Models, Hypergeometric Systems and Dessins d'Enfants. It demon- 
strates that the setting in which Gelfand, Kapranov and Zelevinsky 
have formulated the theory of hypergeometric systems, provides also 
a natural setting for dimer models. The Fast Inverse Algorithm of 
[Hj and the untwisting procedure of [Ij are recasted in this more nat- 
ural setting and then immediately produce from the quiver data the 
Kasteleyn matrix for dimer models, which is best viewed as the bi- 
adjacency matrix for the untwisted model. Some perfect matchings 
in the dimer models are direct reformulations of the triangulations in 
GKZ theory and the rule which maps triangulations to the vertices of 
the secondary polygon extends to a rule for mapping perfect match- 
ings to lattice points in the secondary polygon. Finally it is observed 
in many examples and then conjectured to hold in general, that the de- 
terminant of the Kasteleyn matrix with suitable weights becomes after 
a simple transformation equal to the principal ^-determinant in GKZ 
theory. Illustrative examples are distributed throughout the text. 



1 Introduction 

In the last decade interesting correspondences were discovered relating Quiver 
Gauge Theories, lattice polygons and Calabi-Yau singularities. The motiva- 
tion and evolution of these ideas in physics are well-documented in many 
articles; e.g. [Il[3ll2[5l[6l[T2l[l3l[Tl[El[16l[I71[T9l[20^ In the present paper 
we want to put some aspects of these correspondences from physics alongside 
the hypergeometric systems in two variables of Gelfand, Kapranov, Zelevin- 
sky [HI El lO] and the dessins d'enfants of Grothendieck et al. [T8| [22| [2i] . 
This reveals intriguing connections between these fields. 
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The beautiful insight of Gelfand, Kapranov and Zelevinsky was that hy- 
pergeometric structures greatly simplify if one introduces extra variables and 
balances this with an appropriate torus action |8l [9l |10l Ull |25] . In order to 
profit from the simplication they developed tools like the secondary fan, sec- 
ondary polytope and principal A- determinant. This paper demonstrates that 
these are also very practical tools for studying quivers and dimer models. 

A similar beautiful insight of simplification by going to higher dimensions 
appeared in De Bruijn's construction [2J of Penrose tilings and developed into 
the well-known projection method in the theory of quasi-crystals; e.g. [23] . 
We apply the same method to construct periodic rhombus tilings of the plane, 
which the physicists call brane tilings and dimer models. 

From a geometric perspective this paper deals with embeddings of quivers 
into compact oriented surfaces without boundary. More specifically, the initial 
combinatorial data for a quiver Q are two finite sets E (arrows) and V (nodes) 
and two maps s,t : E ^ V (source and target). Embedding Q into a compact 
oriented surface without boundary Ai means that V becomes a subset of Ai 
and an arrow e & E becomes a path Pf, in M. from the point s(e) to the point 
t[e). It is required that the boundary of every face of (A1, Q) - i.e. connected 
component of \ IJeeE^e - is formed by a sequence of paths {jpe^, . . . ,Pe„) 
with n Pe,+i = t{ei) = s(ej+i) for 1 < i < n; e^+i = ei. It is evident 
from this requirement that the face boundaries receive an orientation from 
Q. When this is combined with the orientation on Ai faces lie either on the 
positive or on the negative side of their boundary. In pictures we mark faces 
which lie on the positive (respectively negative) side of their boundary with 
• (resp. o). It is obvious that adjacent faces get different colors. The dual 
graph of Q w.r.t. Ai is thus a bi-partite graph r*°, embedded in Ai: the 
nodes of r'° correspond with the faces of (A^,Q); the edges of r*° connect 
nodes coming from adjacent faces and correspond bijectively with the edges of 
Q. A pair {A4, Q) consisting of a compact oriented surface without boundary 
and an oriented graph embedded in it is called a dessin d'enfants or just 
dessin, in one of the various equivalent definitions of "dessin (d'enfants)"; see 
[221 [ISl El] • Other definitions refer to {Ai, r*°) as dessin d'enfants. In case Ai 
has genus 1 one often calls r*° and its lifting to the plane (i.e. the universal 
covering of the torus A4) a dimer model or brane tiling; see for instance 
[H El Il2l [131 [in [IS] ■ Yet another presentation of the same structure gives the 
dessin as the triangulation of Ai with vertex set the union of the nodes of Q 
and r*°; the triangles are given by the triplets of vertices consisting of two 
nodes of r*° connected by an edge e* of r*° and one node of Q incident to 
the edge e of Q which is dual to e*. Figure [T] shows these three get-ups of the 
dessin d'enfants of the P^-quiver (case Bi in Figure [3]) 

From an algebraic perspective this paper deals with superpotentials for 
quivers. There is a simple equivalence between the geometric and algebraic 



2 



Figure 1: Three versions of the dessin d ' enf ants for F"^ (= case Bi in Fig. 
The surface is obtained by identifying opposite sides of the hexagon. 

perspectives: the superpotential is a convenient way of writing the hst of ori- 
ented boundaries (pen ■ ■ ■ ,Pe„) of the faces of (A^, Q) with ±-signs indicating 
how the orientation matches with that of A4. It can also be read as the 
instruction for building Ad by glueing polygons. For the dessin d'enfants in 
Figured] the superpotential (for some numbering of the arrows of Q) is 

X1X2X3 + X^X^Xq + XyXsXg — X2XsX^ — X^XqXq — X^X\X-j . 

Another algebraic perspective, equivalent to the previous one, is given by 
what we want to call the bi-adjacency matrix of the dessin d'enfants {Ai, Q) 
with weight w. It is defined as follows. The weight is a function w : E ^ C 
The dessins in the present paper have as many black faces as white faces 
and every edge e E E lies in the boundary of a unique black face, denoted 
b(e), and a unique white face, denoted w(e). The bi-adjacency matrix of 
[Ai, Q) with weight w is a square matrix with rows corresponding with 
the black faces, columns corresponding with the white faces and entries in 
the polynomial ring C[mi, | f G which has one variable for every node v of 
the quiver Q: the {h,w)-entry ofW^ is 

^b,w = Yl Us(e)Ut(e) ■ (1) 

e£E: b(e)=b, w(e)=w 

The bi-adjacency matrix is in fact the Kasteleyn matrix of a twist of the dimer 
model (7W,r*°); see Section [HI For the dessin in Figure [1] the bi-adjacency 
matrix (for some numbering of the nodes, edges and faces) is 

WiUiU'i ZU2U1U2 
W4U2U3 W^UiU-^ 
ZUjUiU2 W^U2U-3 WqUiU^ 

In Section [2] we describe how the quivers for which we can solve the em- 
bedding problem, are associated with certain rank 2 subgroups L of Z^; here 
N is the number of nodes of the quiver. Such subgroups L C are the 
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foundations for the theory of Gelfand, Kapranov and Zelevinsky. In Section 
[3] we describe the secondary fan and the secondary polygon associated with 
L C Z^. These important structures in GKZ-theory surprisingly turn out to 
be also quite relevant for the dimer models. We show for instance in Theorem 
19.31 that the secondary polygon is the Newton polygon of the determinant of 
the bi-adjacency matrix with non-zero weights. In Section [10] we report the 
observation, based on examples, that for the critical weight 

crit : ^ ^ Z>o , crit(e) = ^{e G E \ s(e') = s(e), t{e') = t{e)} (2) 

the determinant of the bi-adjacency matrix K'^'''*, after a simple transforma- 
tion, becomes equal to the principal A-determinant of Gelfand, Kapranov, 
Zelevinsky [11] . In Sections H] and [5] we recall some topics in the theory of 
GKZ-hypergeometric systems of differential equations. 

Section E] presents the algorithm to solve the quiver embedding problem. 
It is a combination of the Fast Inverse Algorithm of |14j and the untwisting 
procedure of [1]. In the cited papers, however, methods are mainly presented 
via visual inspection of pictures in some concrete examples. So an extrapo- 
lation to general situations was needed. The first impression that the Fast 
Inverse Algorithm of [Hj is more or less De Bruijn's construction [2] of Pen- 
rose tilings did not quite yield the sought after embedding of the quiver; the 
untwisting procedure of ^ is also needed. The algorithm became a smoothly 
operating algebraic-combinatorial tool by consistently working from the phi- 
losophy that things look simpler from a higher dimensional viewpoint (with 
group action). For doing the computer experiments behind this paper I imple- 
mented the algorithm in MATLAB. The reader can find in Section [6] sufficient 
details for making a computer version of the algorithm. 

Sections [71-[T0l demonstrate that consistently working from the higher di- 
mensional viewpoint leads to new insights in the dimer technology tools per- 
fect matchings, Kasteleyn matrix and its determinant and shows their close 
relation with the GKZ tools secondary polytope and principal A-determinant. 

Besides problems such as proving that the algorithm of Section [6] yields 
at least one superpotential for every quiver which satisfies the conditions in 
Theorem 12.101 or proving Conjecture 110.51 our work raises some interesting 
questions like: 

Ql. The Calabi-Yau singularities in the background of this work are 
constructed from toric diagrams, which here are interpreted as secondary 
polygons. In [19] §6 the singularity is obtained from the toric diagram by 
a symplectic quotient construction. In [12] §2 the singularity is given as a 
toric variety constructed from a fan with one maximal cone, namely the cone 
over the toric diagram (polygon). For this the polygon is put into the plane 
of points with first coordinate 1 in M'^. In the work of Gelfand, Kapranov 
and Zelevinsky the secondary polygon is put into a 2-dimensional plane in 
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M , which does not pass through 0; here is the number of nodes of the 
quiver. One can perform the standard toric variety construction with the 
3-dimensional cone over the secondary polygon in M^. Does this toric variety 
give the appriopriate view in the GKZ philosophy (higher dimension compen- 
sated by group action) on the singular Calabi-Yau 3-space? Since this toric 
variety is a natural domain for GKZ hypergeometric functions (see 15.21) one 
may wonder: Do hypergeometric functions provide new useful tools for in- 
vestigating the geometry of 3-dimensional Calabi-Yau singularities? Possibly 
positive indications are the fact that the dimension of the solution space of 
the GKZ system of hypergeometric differential equations equals the size of 
the bi-adjacency matrix K'^'''* (see Theorem 15.31 and Corollary 17.31) and the 
observation (Conjecture I1U.5I) that the determinant of W^^^'^ equals, up to a 
simple transformation, the principal ^-determinant, which describes the sin- 
gularities of the GKZ system. 

Q2. The surface Ai with embedded in it the pair of dual graphs Q and 
r*° has been constructed in a purely combinatorial topological way from the 
superpotential. It is a general fact (see [TH] §1-2) that there is then a con- 
tinuous map / from Ai to the 2-sphere S*^ which is a ramified covering with 
exactly three ramification points such that the set of vertices of Q, the set of 
black vertices of r*° and the set of white vertices of r*° are the fibers over the 
three ramification points. A highlight in the theory of dessins d'enfants is Be- 
lyi's Theorem (see [TS] Theorem 2.1.1). It states that in the above situation 
the surface A4 admits a model Aix over a number field K such that / be- 
comes a morphism Aix P^- of algebraic curves over K which is unramified 
outside {0, 1, oo}. The labeling of the ramification points can be taken such 
that /^^(oo), f^^{0) and f^^{l) are the sets of white and black vertices of 
r*° and the vertices of Q, respectively. The bipartite graph r*° in Ai is then 
the inverse image of the negative real axis [— oo,0] in the Riemann sphere 
Pj. = and the quiver Q is the inverse image of the positively oriented 
unit circle {2;gC||2;| = 1}. It is usually difficult to find explicit algebraic 
equations for Aix and the Belyi function f. 

On the other hand, the authors of |4j write in footnote ^: .... we have pro- 
duced a dimer model that is defined on its own spectral curve det = 0. 
Unfortunately their arguments are not (yet) sufficiently refined to yield the 
weight UJ that is to be used in this equation. 

It seems an interesting challenge to tackle the two problems simultaneously 
and look for a weight zu and a Belyi function / on the algebraic curve with 
equation det = that realize {Ai, Q, r*°) as described above. 
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2 Rank 2 subgroups of Z and Quivers. 



2.1. The most basic object in tliis paper is a rank 2 subgroup L C Z 
wliicli is not contained in any of the standard coordinate hyperplanes of JJ^ 
and is perpendicular to the vector (1, . . . , 1), with aU coordinates 1. So, the 
elements of L are vectors (£i, . . . , ^at) G with £i + . . . + f^v ~ and for 
every i e {1, . . . , A^} there is an (£i, . . . , £jv) G L such that ii ^ 0. 

2.2. Notation. Throughout this paper ei, . . . ,eiv is the standard basis of 
and J denotes the matrix . 

2.3. Definition/construction. Taking the second exterior powers one 
finds an inclusion /\^ L "-^ /\^ Z^. The group /\^ L is a free Z-module of rank 
1. After fixing the orientation on L, i.e. choosing one of the two possible 
isomorphisms /\^L ~ Z, one gets an inclusion Z ^ /\^ Z^. The coordinates 
of 1 e Z with respect to the standard basis {ei A Gj}i<i<j<N of A 

are 

called the Pliicker coordinates o/L C Z^. 

Dual to the above inclusion one finds the map /\^ Z^ Z, i.e. an anti- 
symmetric bilinear form on Z^. In view of its relation with the Pliicker 
coordinates we call this the Pliicker form of L. 

To make this look more explicit we take a basis for L compatible with the 
chosen orientation. One can represent the two basis vectors as the rows of 
a 2 X A^-matrix B. Let bi, . . . , b^r be the columns of this matrix. Then the 
Pliicker coordinates are (det(bj, bj))-^^^^^.^^ and the (anti-symmetric) matrix 
for the Pliicker form is (det(bj, bj))i<ij<iv- The Pliicker coordinates and the 
Pliicker form do, of course, not change if one takes another basis for L with 
the same orientation. Note that the matrix for the Pliicker form has rank 2: 

rank ((det(bi, bj))i<ij<iv) = 2. (3) 



2.4. Definition. A quiver Q is a finite directed graph. It can be given 
combinatorially by two finite sets E and V and two maps s,t : E ^ V; in 
short hand notation Q = {E,V,s,t). Each element of F is a node of the 
graph and an element e e £J is an arrow from node s(e) to node t{e). The 
adjacency matrix of the quiver is the matrix {qij)ij,^v with entry q^j equal to 
the number of arrows from node i to node j. If the quiver has no directed 
loops of length < 2 (i.e. qijqji = for all i, j E V), it is faithfully described by 
the anti- symmetrized adjacency matrix Q = {qij — qji)ij^v If the (i,j)-entry 
of Q is positive it equals the number of arrows from node i to node j. If it is 
negative it is the number of arrows from j to i. 
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2.5. Definition. We call the quiver without directed loops of length < 2 
with anti-symmetrized adjacency matrix equal to the matrix of the Pliicker 
form of L the Pliicker quiver of L. 

2.6. Examples of groups L C and their Pliicker quivers are shown in 
Figures [3] and IH 

2.7. Some easily visible geometric properties of a quiver without directed 
loops of length < 2 are equivalent to easily seen algebraic properties of its 
anti-symmetrized adjacency matrix Q. For instance, the graph has no isolated 
nodes if and only if no column of Q is zero. Also, in every node of the 
quiver the number of incoming arrows equals the number of outgoing arrows 
precisely if in each row of Q the sum of the entries is 0. For Pliicker quivers 
these properties correspond to the conditions that L is not contained in any 
of the standard coordinate hyperplanes of and is perpendicular to the 
vector (1, . . . , 1). One property of a Pliicker quiver which one can not easily 
read off from the graph is rank Q = 2. 

In the remainder of this section we show that these properties characterize 
Pliicker quivers of rank 2 subgroups L C as in 12.11 

2.8. Proposition. An N x N -matrix C with entries in Z is anti- symmetric 
and rank C = 2 if and only if there is a 2 x N -matrix B with entries in Z 
such that rank B = 2 and C = B^ J B. 

Proof. The "if'-statement is trivial. So let us consider the "only if" and 
assume that C is anti-symmetric and rankC = 2. Let d denote the greatest 
common divisor of the entries of C and C = -.C. Choose an x 2-matrix 

d 

D whose columns form a Z-basis for the column space of C. The equality of 
column spaces of C" and D means that there are an x 2-matrix E and a 
2 X A^-matrix F, both with entries in Z and of rank 2, such that D = C'E 
and C = DF. Then C = DF = C'EF = -C'^EF = -F^D^EF . So 
D^E is an anti-symmetric 2 x 2-matrix; say D*E = /J with f E 1j. In fact 
/ = ±1 since it divides all entries of C. If / = 1 we replace F by ^F. We 
then always have C = F^SF. Let G be any 2 x 2-matrix with entries in Z 
and det G = d. Let B = GF. Then C = B^SB as wanted. ■ 

2.9. Multiplying in Proposition 12.81 the matrix B from the left by a matrix 
from (Z) does not change the matrix B^SB. So it is more natural to 
interpret C as the matrix of the Pliicker form of the Z-row space of B. Note 
that it follows from the proof of Proposition 12.81 that this space is uniquely 
determined if the greatest common divisor of the entries of C is 1. If on the 
other hand the greatest common divisor of the entries of C is d > 1 this 
space depends on the additional choice of a 2 x 2-matrix with entries in Z 
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and determinant d. 

Let us summarize the above discussion: 

2.10. Theorem. For every L C as in \2.1\ one has its Pliicker quiver. 
This quiver has no isolated nodes or directed loops of length < 2 and in every 
node the number of incoming arrows equals the number of outgoing arrows. 
Moreover the rank of its anti- symmetrized adjacency matrix is 2. Conversely, 
every quiver with these properties is the Pliicker quiver of some L C as in 
\2.1\ The correspondence between such quivers and such L C Z^ is one-to-one 
for those quivers for which the greatest common divisor of the entries of the 
anti- symmetrized adjacency matrix is 1 and those L C Z^ for which Z^L 
has no torsion. ■ 



3 The secondary fan and polygon of L C Z . 

3.1. For L C Z^ as in 12. II we now define the secondary fan and the secondary 
polytope. In 14. 121 we will compare this with the original definitions by Gelfand, 
Kapranov and Zelevinsky. The term "secondary" refers to the fact that in 
their theory of hypergeometric systems another polytope appears first, which 
is therefore called the primary polytope. Nonetheless both the secondary fan 
and the secondary polytope can most conveniently and directly be described 
using the lattice L. In case the rank of L is 2, the constructions become 
particularly simple. 

3.2. Definition. Let L C Z^ be as inO Let := Hom(L,]R) denote the 
real dual space of L. Let ei, . . . , e^r be the standard basis of Z^. Let bj G 

be the image of e^ under the map dual to the inclusion L ^ Z^. 

Here and henceforth we identify with the real dual space of Z^ by means 
of the standard dot product. 

By definition, the secondary fan ofh is the following collection of cones in 
L^: the 0-dimensional cone {0}, the 1-dimensional cones M>obj (z = 1, . . . , A^) 
and the 2-dimensional cones which are the closures of the connected compo- 
nents of L]^ \ [jf^^ M>obi. 

3.3. Definition. With a 2-dimensional cone C in the secondary fan one 
associates the set 

Lc ■■= {{^, J}C{1,...,N}\CC (M>ob, + R>ob,) } (4) 
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of 2-element subsets of {1, . . . , A^} and the vector 

i>c ■= I det(bi, hj) \ (ci + e^-) . 

The secondary polytope S(L) of L C is then defined as 

S(L) := convex hulld-i/^c I C 2-dim cone of secondary fan}) . 



(5) 



3.4. Remark. In §4.121 we give an interpretation of the above Lc as a 
triangulation of the primary polytope. In Equation (l27|l we associate with 
Lc a perfect matching in a bipartite graph. 

3.5. Example. Figure [2] shows the secondary fan with the lists Lc, the 
secondary polytope S(L) with coordinates for the vertices and the primary 
polytope for L = Z(0,l,l,-2) © Z(-l, 0,2,-1) C Z\ This is case B2 of 
Figure [31 see also Example 14.61 and §4.121 
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Figure 2: The secondary fan (right) and the primary polytope (top left) and the 
secondary polytope (bottom left) for h = Z(0, 1, 1, — 2)©Z(— 1, 0, 2, — 1) C Z^. 



3.6. The geometric formulation of the construction of the secondary fan 
and polytope is quite attractive. Nonetheless in practical algorithms every- 
thing can be most easily obtained from the matrix of the Pliicker form of 
L C Z^. Indeed, a 2-dimensional cone in the secondary fan is bounded by 
two half-lines M>obj and M>obj with the property that det(bj,bj) > and 
det(bj, bfc) det(bj, b^) > for A; = 1, . . . , A^. Thus to find the secondary fan 
one just needs to find all pairs b^, hj with these properties. 
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For a cone C bounded by ]R>obj and M>obj the list Lc is then 

Lc = { {fc, /} I det(bfc, hi) > and det(bj, b/) > } . 

Now consider two adjacent 2- dimensional cones in the secondary fan, say 
C and C. Let bj^ , . . . , bj^ be those vectors from the set {bi, . . . , b^v} that lie 
on the half-line C H C. Without loss of generality we may assume C lies to 
the right of bj^ and C to the left. Then 

LcMLcnLc) = {{kj}\l e {ii,...,is} , det{hi.hk) <0} , 
Lc'MLcnLc) = {{kj}\l e {ii,...,is} , det{hi,hk) >0} . 

From this we see, using det(bj, b^) = for all i, that 

s N 

i^c - ipc = ^^det{hi^,hk)ek . (6) 

r=l k=l 

In other words, i/jc' — i^c is the sum of rows ii, . . . ,is of the matrix of the 
Pliicker form. 

3.7. To make this even more explicit and simple looking we take a basis for 
L compatible with the chosen orientation. One can represent the two basis 
vectors as the rows of a 2 x A^-matrix B. Let bi, . . . , bjy be the columns of 
this matrix. The elements of L should now be written as row vectors with 
two components and the embedding L ^ TL^ is given by v X]fcLi(v-bfc)efc. 
Since det(bj,bfc) = b*Jbfc, we can reformulate Equation ([6]) as 

s 

■0(7/ — '^c = image of b*^ J under the embedding L . (7) 

r=l 

In order to obtain the simplest formulation for the construction we order 
the vectors bi, . . . , bAr so that the points pfe = Y^i=\ b* for A; = 1, . . . , lie 
ordered counterclockwise on the boundary of the polygon 

A = convex hulljpi, Pat} . (8) 

Then the secondary polygon S(L) is obtained by first rotating A clockwise 
over 90° , next embedding it along with L into and finally translating it 
over the vector ipc, where C is the cone in the secondary fan with left hand 
boundary R>obi. 

Note that while the secondary polytope S(L) depends only on the embed- 
ding L ^ Z^, the polygon A usually changes when one puts another vector 
bj in first position by a cyclic permutation or when one multiplies the vectors 
bi, . . . , b^v by a matrix from S'/2(Z). 
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3.8. For a converse to the above construction we start from a convex polygon 
A in with vertices in 1? and non-empty interior. Let denote its 
boundary. Next we choose a collection of points pi, . . . , pat in fl 1? which 
includes all vertices of A. We number these points so that pj and pi+i are 
consecutive points in the counter-clockwise orientation of 9A. Thinking of 
the elements of 1? as row vectors we define column vectors bi = p* and 
bj = p* — \>\_x for i = 2, . . . ,N . Finally we define L to be the Z-row space of 
the 2 X A^-matrix B with columns bi, . . . , bAr. 

Thus every convex polygon A in with vertices in 1? and non-empty 
interior can he viewed as the secondary polygon of some rank 2 subgroup 
L C satisfying the conditions in \2.1\ 

Note however that in general there are several possible choices for the 
points pi, . . . , ptv in (9A fl 1?. The minimal choice takes only the vertices of 
A, while the maximal choice takes all points of dA fl 1?. 

3.9. Theorem. The (Euclidean) area of the polygon A in ^ is 

area A = \ Ei<i<j<7v I det(bi, b^)] - E{i,j}gLc I det(bi, b^)] (9) 
for every 2-dimensional cone C in the secondary fan. 

Proof. After applying, if necessary, a cyclic permutation to bi, . . . ,bjv we 
may assume, without loss of generality, that C = ]R>obi + ]R>ob7v. 

From the triangulation A by the diagonals between the vertex pjv = 
and the other vertices of A one sees that area A = | ^-^^^^^.^^^ det(bj, bj) 
where M is such that b^ ^ ]R>ob7v for 1 < j < M and hj G ]R>ob7v for 
M + 1 < j < N. The fact that Y.f=i = 0^ implies 

X]i<j<j\/<j<Ar "^^^bj, bj) = det(^-|^<^<j^^ bj, bj) = 0. 

Thus we see 

areaA = | Ei<i<j<7v det(bi, b^) . 

Next note that for 1 < z < j < the inequality det(bj,bj) < holds if 
and only if C C M>ob,; + M>obj, if and only if {i,j} G Lq- Equality ([9]) now 
follows immediately. ■ 

3.10. Corollary. The number of lattice points in the interior of the polygon 
A in ^ IS 

-2 {^i<i<j<N I det(b„ b,)| - iv) - E{^,}eLc I det(b„ b,)| + 1 . 

Proof. This follows from Theorem 13.91 in combination with Pick's Formula 
([H p.ll3): 2areaA = 2jj (Z^ n interior A) + [t (Z^ n boundary A) - 2 ■ 
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Figure 3: Quivers from Figures 10, 11, 4, 12 and corresponding polygons. 
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3.11. Example. Figures [3] and H] show some examples of the above relation 
between quivers, rank 2 subgroups L of TL^ as in 12.11 and polygons. The 
quivers are taken from |3]. We name our examples neutrally as . . . , -Bio, 
but also mention the names given to these quiver models in [3]. The latter 
names refer to del Pezzo surfaces dP^, i.e. the projective plane blown 
up in k points, some of which, on the present occasion, may happen to be 
"infinitely near" and therefore require repeated blow-ups. With the latter 
somewhat liberal use of the name del Pezzo surface, dP^ matches well with 
the toric geometry of the fan whose 1-dimensional rays are the half-lines from 
the interior lattice point o through a lattice point • at the boundary. Note 
however that this fan is in general not the same as the secondary fan and 
that therefore the singularity associated with the polygon (cf. [I2l [19]) is 
not the singularity obtained by contracting to a point the zero-section of the 
canonical bundle of the del Pezzo surface. 

The polygons in Figures [3] and HI however, are in general different from, 
the polygons associated with the same quiver in [3] Figure 8. We will clarify 
this issue in lOl [6J[5l lOGl 

It is remarkable that the polygons in Figures [3] and H] are exactly the 
lattice polygons with one interior lattice point and < 6 lattice points on 
the boundary. There are exactly 16 lattice polygons with one interior lattice 
point. The remaining 6 are shown in Figure [5l A nice review of these so-called 
reflexive polygons can be found in [21] . 

In cases and the greatest common divisor d of the entries in 

the quiver's adjacency matrix is > 1 and the polygon in these cases depends 
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on an integer matrix with determinant d. We have chosen this matrix such 
that the polygon is a reflexive polygon. 

4 GKZ hypergeometric systems 

4.1. In the late 1980's Gelfand, Kapranov and Zelevinsky discovered fas- 
cinating generalizations of the classical hypergeometric structures of Euler, 
Gauss, Appell, Lauricella, Horn P, El UHl ttU The main ingredient for 
these new hypergeometric structures is a finite sequence A = (ai, . . . , a^v) of 
vectors in Z^"*"^ which generates Z'^"^^ as an abelian group and for which there 
exists a group homomorphism h : Z'^"'"^ — >■ Z such that /i(aj) = 1 for all i. 
The latter condition means that A lies in a A;-dimensional affine hyperplane in 
^fc+i_ pigure shows A (each black dot represents one vector) sitting in this 
hyperplane for some classical hypergeometric structures. Gelfand, Kapranov 
and Zelevinsky called these new structures A- hypergeometric systems, but 
nowadays most authors call them GKZ hypergeometric systems. 

Let L denote the lattice (= free abelian group) of linear relations in A: 

L := . . . , £^) G Z^ I £iai + . . . + ij^Sij^ = 0} . (10) 

It follows from the above construction and assumptions that the quotient 
group is Z'^^^, and a fortiori, it is torsion free. 

In this paper we consider only the case when the rank of L is 2; i.e. 
k + 1 = N — 2. Since we assumed the existence of a linear map h : Z^~^ Z 
such that h{a.i) = 1 for all i, L lies in the kernel of the map s defined by: 

s:R^^R, s{zi,...,Zn) = Zi + ... + zn. (H) 

The following lemma shows what the other condition for L in 12.11 means 
in terms of A. 

4.2. Proposition. In the situation of \4-l\ the following statements are equiv- 
alent: 

1. h is not contained in any of the standard coordinate hyperplanes o/Z^. 

2. No vector in the sequence A is linearly independent of the other vectors. 

3. For every i the lattice of linear relations between the vectors of A \ {aj} 
has rank 1. 

Proof. Assume L is contained in the i-th coordinate hyperplane. Then 
the i-th component of every element of L is and hence a^ occurs in no 
linear relation for the vectors in A. In other words, L is the lattice of linear 
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relations between the vectors of ^ \ {aj}. Since L has rank 2, we conclude 
that statement 3 imphes statement 1. 

Conversely, assume that the lattice Lj of linear relations between the vec- 
tors of ^ \ {aj} has rank > 1. Since Lj is contained in L and the latter has 
rank 2, we see that Lj ® M = L CS) M. Therefore, since Lj is contained in the 
i-th coordinate hyperplane, so is L. Hence statement 1 implies statement 3. 

It is obvious that statements 1 and 2 are equivalent. ■ 

4.3. Definition. A subsequence A' of A is said to be minimally dependent if 
the vectors in A' are linearly dependent and the vectors in every subsequence 
A" of A' with A" 7^ A' are linearly independent. A linear dependence relation 
J2j dj^j = is a minimal linear dependence relation in A if the subsequence 
(a^ I aj 7^ 0) is minimally dependent. 

4.4. Proposition. 

1. Let B be a 2 X N -matrix such that its rows are a Z-basis for L. Let 
bi, . . . , b^r be the columns of this matrix. Then one has for every i 

N 

^det(bj,bj)aj- = 0. (12) 



2. Assume the conditions in Proposition \4 . 2\ are satisfied. Then every lin 



ear relation f[T2|) is minimal and every minimal linear dependence rela- 
tion in A is a non-zero scalar multiple of some relation f|T2l) . 

Proof. For statement 1, let ^4 be the matrix with columns ai, . . . , aAr. Then 
Equation flTU]) can be rewritten as BA^ = 0. This implies B^SBA^ = 0, which 
is just a compressed form of the relations we wanted to prove. 

For statement 2 note that condition 3 in Proposition U]2] implies that every 
linear relation (fT2|) is minimal. Conversely, consider a minimal linear relation 

ajBLj = 0. Then the subsequence (a^- 1 aj ^ 0) has at most A^ — 1 elements. 
The given minimal relation is therefore a non-zero scalar multiple of at least 
one relation (1121) . ■ 

4.5. Gelfand, Kapranov and Zelevinsky [8],|9l|T0] associate with a set A as in 
14. H and a vector c G C'^"^^ the following system of partial differential equations 
for functions $ of variables Mi, . . . , tt^y: 

• for every (£i, . . . , Ifq) G L one differential equation 

n(^)"'*- n(0'*. (13) 
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the system of + 1 differential equations 



9^ 

aiui- h... + a7vM7v- — = c$. (14) 

OUi OUn 



4.6. Example. The sequence A 
Figure [3] can be taken to be 



(ai, a2, as, a4) G for Example B2 in 





' 1 " 




' 1 " 




' 1 " 




" 1 " 


ai = 





, ag = 


1 


, a4 = 


2 


, a2 = 


3 



It is a classical result of K. Mayr that the roots of the 1- variable cubic 
polynomial P{x) = ui + u^x + 1/40;^ + U2X^ as functions of the coefficients 
Ui, . . . ,U4 satisfy the GKZ system of differential equations for this A and 

^ ' ; see e.g. [25] §2.1. 



with c 



4.7. Example. An important example in the toric geometry constructions 
of Sasaki- Einstein manifolds is known under the name L"'''''^; see e.g. [6]. 
Here a, b, c are integers with c < b and < a < 6 The polygon for this 
example, displayed in [6j Figure 2, is a quadrangle with vertices (0,0), (1,0), 
{ak, b), {—al, c) where k and / are integers such that ck + bl = 1. The method 
explained in 13.81 yields L = Z(l, ak — 1, —al — ak, al) ©Z(0, b,c — b, —c). From 
this we see that we can take A = (ai, a2, as, a4) G with 



ai 



1 




' 1 " 




" 1 " 




" 1 " 




, a2 = 




, as = 





, a4 = 


b 


c — a 




c 





As in Example 14 . 61 this means that the corresponding GKZ system deals with 
the roots of the "four-nomial" Uix'^~°' + U2x'^ + M3 + u^^x'^ as functions of the 
coefficients ui,U2,u-i, u^. 

According to [6J §3 the examples Y^''' of [H US] are special cases of the 
above: Y^''^ = L^-i^p+i^p, So, these correspond to the "four-nomials" 
Uix'^ + U2x''^ + U3 + UiX^^'^. Thus Example 14.61 is in fact Y'^'^ = L^'^'^. 

4.8. Example. From the pictures of ^ shown in Figure [6] one easily sees that 
for the Gauss system the lattice L is generated by the vector (1,1,-1,-1) 
and thus has rank 1 and is not of interest here. For Appell's Fi and F4 the 
lattice L has rank 2 and the corresponding S-matrices are 



Fi : 



1 -1 
1 



0-110 
-1-10 1 



Fa 



-1-1110 
-1-10 11 
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Figure 6: The sets A for Gauss's hypergeometric functions (left), Appell's Fi 
(middle) and Appell's F4 (right). 

Thus Appell's Fi corresponds to Bj in Figure HI For Appell's F4 the 
quiver and the polygon are 



3 




In [8] §3.2 Gelfand, Kapranov and Zelevinsky discuss how the 14 complete 
hypergeometric Horn series in two variables fit their theory. It is an amusing 
simple exercise to now find the corresponding quivers and polygons. It turns 
out that, with the exception of F4, all Horn series give a quiver listed in 
Example 13. lit 

{G3}^B2, {H,}^B,, {Gi,H3,He}^B,, {F,, G2} ^ B, , 
{H^}^Bs, {H,,Hr}^B,, {F^, F3, H^} ^ Bn ■ 



4.9. In 14. II we started from the sequence of vectors A and then defined L via 
Equation (ITU]) . Let us reverse the procedure and start with a rank 2 subgroup 
L C as in 12.11 We can now define 

A = (ai,...,a7v) with (15) 
aj = the class mod L in the quotient group ■ 

If the quotient group is torsion free, it is isomorphic to Z^~^ and the 

sequence of vectors A in f|T5l) satisfies the requirements of I4.1[ Thus / HOj) and 
( [73]) give equivalences, inverse to each other, between the data A and L. 
In combination with Theorem 12.101 this gives 

4.10. Theorem. There is an equivalence of data between on the one hand 



sequences A as in\4.1\ which satisfy also the conditions in Proposition \4 .21 and 
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on the other hand quivers without isolated nodes or directed loops of length 
< 2, such that in every node the number of incoming arrows equals the number 
of outgoing arrows and such that for its anti- symmetrized adjacency matrix 
the rank is 2 and the greatest common divisor of the entries is 1. ■ 

4.11. Remark. We refer to 15.51 for a discussion of how to extend the 
definition of GKZ systems so as to allow for torsion in the group ^^L- 

4.12. In order to have an efficient method to construct bases for the solution 
space of the hypergeometric system (fT3|) - (fT^ Gelfand, Kapranov and Zelevin- 
sky developed the theory of the secondary fan and the secondary polytope. 
The term "secondary" refers to the habit of considering the set A as primary 
data and to call the convex hull of A the primary polytope. 

In case the lattice L of relations in A has rank 2, the constructions become 
particularly simple. The secondary fan is the one in Definition 13.21 With a 2- 
dimensional cone C in the secondary fan one associates a set Lc of 2-element 
subsets of {1, . . . , A^} as in Equation (jlj). One can interpret this Lc is as a 
triangulation of the primary polytope convex hull(^) as follows. Recall that 
A = (ai, . . . , sln) and assume for convenience of imagining pictures that all 
aj's are different. For {i,j} G Lc set T^ij} '■= convex hull({afc \ k i,j }). 
This T{j j} is an (A^ — 3)-dimensional simplex and the simplices T^ij} with 
{i,j} G Lc together constitute a triangulation the primary polytope. 

In [TT] p. 220 the secondary polytope is constructed (and defined) in terms 
of these triangulations as follows. For a 2-element subset {i,j} C {1, . . . , N} 
with i < j let Aij denote the (A^ — 2) x (A^ — 2)-matrix with columns a^ for 
k ^ i,j in the natural order of increasing indices. To a 2-dimensional cone C 
in the secondary fan one then associates the vector 

(pc := ^ IdetAjl^efc. (16) 

Then in [11] p. 220 the secondary polytope ^{A) is defined as 

S(^) := convex hulld^j^; | C 2-dim cone of secondary fan}) . (17) 

As I det Aij\ is (A^ — 3)! times the Euclidean volume of the simplex T^ij}, the 
number vol^ := J2{i.j}eLc I every C equal to (A^ — 3)! times the 

Euclidean volume of the primary polytope convex hull(^) and the point 

(Ar-2)voU S I det A.I 

^ ^ ^ {i,j}eLc k^ij 

coincides for every C with the barycentre of the primary polytope. The 
secondary polytope therefore lies in a 2-dimensional plane parallel to Lr. 
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4.13. Let us compare ^{A) with the polytope E(L) of Definition 13.31 It 
is a well-known and easy to prove fact that | det Aij\ = \ det(bj, hj)\ for all 
{i,j} G Lc (see e.g. |25] Eq. (62)). Thus Equations f|T6l) and ([5]) yield 

N 

= -i'c + voU ■ ^efc , (18) 

k=l 

with 

voU = Yl |detA,,|= |det(b,,b,)|. (19) 

This means that the two secondary polytopes S(^) and E(L) are related by 
a point symmetry with centre |vol^ ■ ^k- 

N 



S(^) = - S(L) + voU ■ XI • 

fc=i 

5 About the solutions to GKZ systems. 



5.1. The variables in GKZ theory (see Section Hj) are the natural coordinates 
on the space C"^ := Maps(^, C) of maps from A to C. The torus T'^+i : = 
Hom(Z*''+^, C*) of group homomorphisms from Z'^^^ to C*, acts naturally on 
C"^ and on the functions on this space: for a G T^^'^^, u G C"^, a E A and 
$ : C-^ ^ C: 

{a ■ u)(a) = a(a)u(a) , ($ ■ a){u) = $(a ■ u) . (21) 

The GKZ hypergeometric functions associated with A and c are defined 
on open domains in C"^. One easily sees that if a function $ on C"^ satisfies 
the differential equations (fT3l) then for every a G T'^''"^ the function $ ■ a also 
satisfies these differential equations. So, the torus T'^^^ acts on the solution 
space of the system of differential equations ( fT3l) . 

On the other hand, if $i and $2 are two functions which satisfy the dif- 
ferential equations IHM with the same c, then their quotient = satisfies 

aiUi- \- . . . + ajyUNT: — = 0. 

OUi OUn 

The latter equation is equivalent to being T'^ "''^-invariant. Thus we find: 
All quotients of pairs of solutions of the system l[T3\)-[l^\) are functions 



on simply connected open subsets of the orbit space C'^/f^+'^. In the case of 
interest in the present paper the dimension of this orbit space isN — k — 1 = 2. 
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For c G Z'^^^ the differential equations fll4p are equivalent with $ trans- 
forming under the action of T''~^^ according to the character given by c: 

$.a = a(c)$. (22) 

In particular for c = all solutions of f|T4l) are T'^ "''^-invariant. 

5.2. The space (C*)-^ := Maps(^, C*) of maps from A to C* is a torus of 
dimension which contains T^"^^ as a subtorus. The action fl^T]) of cr G T'^^^ 
on C"^ obviously restricts to the action of T'^"'"^ on (C*)"^ by multiplication. 
The quotient space (C*)-^T'=+i, which is a subspace of C-^f^+i, is the torus 
Hom(L,C*) of group homomorphisms from L to C*. The toric variety as- 
sociated with the secondary fan (cf. Definition 13.21) gives a compactification 
of the torus Hom(L,C*); see [7] for the general theory of toric varieties. It 
is an essential part of the GKZ philosophy that quotients of hypergeomet- 
ric functions should be viewed as being defined on open subsets of this toric 
variety. 

5.3. Theorem, (cf. [8] Theorems 2 and 5, [8]', [26] Prop. 13.5 ) Let 

NA = {xiSLi + . . . + xnsln e ]R^+^ I Vxi G Z>o} , 
ZA = {xiai + . . . + xn^ln e R^^'^ I Vxi G Z} , 
pos(^) = {siai + . . . + xnSln e R''^'^ I G R>o} . 

Assume NA = Z^npos(^), then the dimension of the space of solutions of 
the system of differential equations I[T^) -[I4^ at a general point o/C'^ equals 
the number vol^ in Equation ( fTPj) . ■ 



5.4. Remark. As shown in the proof of Prop. 13.15 the condition 
NA = Z^npos(^) in the above theorem is satisfied if A admits a unimodular 
triangulation. The latter condition is equivalent to: there is a cone C in the 
secondary fan such that | det(bi,bj)| = 1 for all {i,j} G Lc (see Definitions 
EJandESD. 

5.5. Remark. From the discussion in 15.11 one may get an idea about ex- 
tending the GKZ system f lT3|) - f|T^ to the situation in which Z^l has a 
non-trivial torsion subgroup (Z^l)^^^^. Equation ([13]) still makes sense m 
this more general situation, but Equation f|T^ must be adapted. 

Recall from f[T^ a definition of A which also works in the torsion case. 
Let := Hom(Z7L,C*). This is a commutative algebraic group of which 
the connected component of the identity is = T'^"'"^ = IIom(Z'^"^^, C*) 
and the group of connected components is the finite abehan group Ga/ Q\ = 
Hom((Z7L)^^^^ ,C*) ~ (Z^l)^^^^. Formula ^ defines an action of on 
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C"^ and on the functions on C"^. It is now clear how to adapt Equation fll4p : 

let Sj denote the projection of in the free part Z^^^ of and replace 
m by 

_ a$ _ a$ 

a-lUi- \- . . . + ajyUNT: — = c<I> 

OUi OUn 

plus the requirement that the solution should transform according to some 
character of the finite abelian group (^^Lj^^j.^- 

5.6. Example. In case i?i in Figure[3]one has L = Z(2, — 1, — 1)©Z(1, 1, — 2) 
and Z')'L = Z © Z/ 3^. As a generator for the torsion subgroup we take g = 
(1, —1, 0) mod L. Then the polynomials $0 = ^U3+uiU2, $1 = ^u^+uiu-s and 
$2 = l^f + 'U2M3 satisfy the differential equations (fT3|) for L, the differential 
equations ([H]) for c = 2, while g ■ = e^''*''/^^^ for r = 0, 1, 2. 



6 From rank 2 subgroups of Z to Dessins. 

We are going to describe a construction which associates with a rank 2 sub- 
group L of Z^, as in 12.11 dessins d'enfants, i.e. bipartite graphs embedded 
in oriented Riemann surfaces. 

The construction is given in [13] §5 and §6 as Fast Inverse Algorithm. In 
op. cit., however, this algorithm is only presented via explicit visual inspec- 
tion of pictures in some concrete examples. In this section we want to present 
a general principle behind the Fast Inverse Algorithm of [H] which uses only 
linear algebra and can be performed by computer. 

The main part of our construction produces rhombus tilings of the plane 
and is the same as N.G. de Bruijn's p] construction of Penrose tilings. His 
work also led to effective methods for making quasi- crystals |23j. Since L is 
defined over Z the construction yields in our situation periodic tilings of the 
plane = M^, not just quasi-periodic ones. We can therefore pass to 
modulo the period lattice and find a tiling of the two-dimensional torus. 

6.1. In this section L is as in 12.11 The quotient group is allowed to have 
torsion. Let denote the real 2-plane in which contains L. We write 
A + Lk for the real 2-plane in obtained by translating over a vector 
A G M^. On the other hand one has in the standard N-grid consisting of 
the hyperp lanes 

'Hi,k := {{zi,...,zn) eR^ \ Zi = k} for i = 1, . . . , AT and A; G Z . 

By intersecting with this standard A^-grid we obtain in the 2-plane A + Lr an 
A-grid of lines: 

:=Hi,fcn(A + LM) for i = 1,. . . , A and A; G Z. (23) 
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Note that this crucially uses the assumption that L is not contained in any 
of the standard coordinate hyperplanes in Z^. 

The grid lines are naturally oriented: the orientation on is such that 
the points in A + Lr which are to the left of have a larger i-th coordinate 
than the points to the right of jC^j.. 

6.2. Definition. We say that A is non-resonant if no three grid lines pass 
through one point. 

6.3. First step. Let L and A be as in 16.11 with A non-resonant as in 16.21 
Consider the map 

F : ^ , F{zi, ...,zn) = { [2iJ , . . . [zn\ ) , 

where [z\ for a real number z denotes the largest integer < z. The map 
F contracts an open A^-cube in the structure on (and part of its 
boundary) onto one of its corners. F commutes with the translation action 
of Z^ on and Z^. Suppressing L and A from the notation we define 

5°:=F(A + Lm). (24) 

The map F is constant on each 2-cell ( = connected component) of the grid 
complement in A + Ljr. Every point of iS° is in fact the image of a unique 
such 2-cell. The distance between two points of is 1 if and only if the 
corresponding 2-cells are separated by exactly one grid line. 

Because A is non-resonant, an intersection point of grid lines is in the 
closure of exactly four 2-cells. If x = C^j^ fl Cj^^ the four points of 5° corre- 
sponding to these cells, i.e. F(x.), F{x.) — e^, -F(x) — e^, -F(x) — e — e^, are the 
vertices of a square Dx. When x runs through the set of all intersection points 
of grid lines, these squares fit together to a connected surface S := IJ^ Dx 
embedded in M^. One can characterize the surface S also as 

S := union of all unit squares in with vertices in the set . (25) 

Here one may define unit square as the convex hull of four points Pi, P2, Pa, Pa 
in with Euclidean distances ||pi — P2IHIIP2 — PsIHIIps — P4IH ||P4~Pi||= 1 
and IIpi— P3|| = ||p2~P4||= V^- A unit square with vertices in Z^ is necessarily 
of the form p + convex hull ( 0, ej, e^, ej -|- e^) for some p G Z^ and some 
We say that a unit square in has type {i,j) if its sides are parallel to the 
vectors e.j or ej. Similarly, a side of a unit square in has type i if it is 
parallel to the vector ej. 

6.4. Proposition. Let = Hom(L,]R) denote the real dual ofh. Identify 
the real dual Hom(Z^, M) o/Z^ with M.^ by means of the standard dot product. 
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Figure 7: converting grid to tiling for B2 from Figure 



Then the linear map which is dual to the inclusion L — * Z^, 

restricts to a homeomorphism S — ^ L^. 

Proof. Let bi, . . . , G denote the images of the standard basis vectors 
of M^. After identifying in the obvious way, the real plane L]r with the real 
2-plane A + Lr in one obtains in Lr an A^-grid in which the grid lines 
are perpendicular to bj. After choosing a basis for L and taking the dual 
basis for we may identify Lr and with M?. The term "perpendicular" 
then means perpendicular with respect to the standard inner product on M^. 

Next draw for each intersection point of grid lines, say x = fl -Cj'm, 
a parallelogram with centre x and sides ebj and ehj. Here the positive real 
number e is so small that the parallelograms obtained from all grid intersection 
points are disjoint; see Figure [7] for an example. It is clear from this figure 
how one can glue the parallelograms of two consecutive intersection points on 
a grid line along their sides perpendicular to the grid line. The result of this 
glueing is then the same as the image (scaled with a factor e) of the surface 
S under the projection L^. ■ 

It is obvious from the construction that the group L acts by translations 
on the point set 5° and on the surface S, preserving the types of the squares 
and their sides. The cell structure on S given by the squares, their sides and 
vertices induces therefore on iS/L a cell structure. The following proposition 
counts the cells on <S/i^. 

6.5. Proposition. Let B be a2x N-matrix such that its rows are a Z-basis 
for L. Let bi, . . . , b^r be the columns of this matrix. Then on S/i^ 

1. the number of2-cells of type is | det(bj, bj)|. 

2. the number of 1-cells of type i is J2f=i I det(bj, bj) | . 

3. the number of 0-cells is ^i<j<,<7v I det(bj,bj)|. 
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Proof. The counts can be made on the dual structure, which is given by the 
intersecting grid hues on Lr. Let Bij denote the 2 x 2-matrix with columns bj 
and bj. Then 1. amounts to counting for Aj, Xj G M the number of elements 
in the set {x = (xi, X2) G | < xi, X2 < 1 , xSjj + (Aj, Xj) G Z^}. Since 
there is a bijection between this set and the coset space ^lTl?Bij., the number 
of elements is | det Bij\. The number of 1-cells of type i is equal to the number 
of intersection points of the i-th grid with the other grids. In view of 1. this 
number is therefore as stated in 2. Knowing the numbers of 0- and 1-cells on 
the 2-torus Lr/L one computes the number of 2-cells from the fact that the 
Euler characteristic of a 2-torus is 0. ■ 

6.6. Second step. For the constructions in 16.11 and 16.31 we may randomly 
choose a non-resonant A G M"^. Once we have constructed the set and 
surface S we modify it by elementary transformations. For such an elementary 
transformation we need a point po G such that there are exactly three 
points pi, P2, Ps in 5° at distance 1 from po and such that the three points 
Pi + P2 — PO) Pi + Ps ~ Po and P2 + Ps — Po also lie in 5°. These seven points 
are the vertices of three squares in S. The elementary transformation now 
replaces po by pg := pi + P2 + Ps — 2po and replaces the three squares by 
the three squares in with vertex configuration {pi, p2, Ps, Pi + P2 — PO) 
Pi+Ps~PO) P2 + Ps~Po; Po} j see FigureO In order to preserve L-periodicity 




Figure 8: An elementary transformation 



we perform this transformation simultaneously at all configurations {po + £, 
Pi + ^, P2 + ^, Ps + ^} with ieh. 

After this elementary transformation we obtain a surface S' in which 
is the union of all unit squares with vertices in the set iS'°, obtained from 
by replacing the points Po + ^ by pi + p2 + Ps — 2po + i for all i Eh. 

The surface S', in turn, can be further modified by elementary trans- 
formations: just replace in the above construction by S'^ and po by an 
appropriate point of S'^. 

6.7. Third step. After starting the construction in 16.31 with a randomly 
chosen non-resonant A G one can create by repeated elementary trans- 
formations many surfaces S each of which is the union of unit squares with 
vertices in an L-invariant subset iS° of and which is mapped homeomor- 
phically onto by the projection L^. Moreover the numbers of cells 
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on each of these surfaces are still the same as in Proposition 16.51 

The number of surfaces one can make in this way is finite and it is possible 
to generate (by computer) a complete list. 

6.8. Remark. The list of surfaces produced in 16.71 is. up to permutation of 
its entries, independent of the choice of A G at the start of the algorithm. 
This can be seen as follows. Fix i, 1 < i < N. Take A G and let 
At = A + te, for t G M. Set ti = min{t G M>o | A* is not non-resonant} and 
t2 = min{t G M>tj | At is not non-resonant}. Then the surface St constructed 
from \t in 16.31 is equal to the surface iSq for < t < ti. Figure [9] shows how 
the grid locally changes as t passes through ti. Comparing Figures [3 and [S] 




Figure 9: Local change in grid at resonance 

one sees that such a local change of the grid corresponds to an elementary 
transformation of the surface Sq. Thus the surfaces St with ti < t < t2 are 
obtained from Sq by a number of elementary transformations. 

Since any two non- resonant A, A' G can be moved to a common value 
by coordinatewise changes as above we conclude that the surfaces S and 
S' obtained in 16.31 from A and A', respectively, are related by elementary 
transformations, and that, hence, the lists of surfaces which the algorithm 
produces from start values A and A' are the same, up to possibly a permutation 
of the entries. 

6.9. Definition. We say that a surface S on this list is perfect if the function 
s : M, s{zi, . . . , zjy) = zi + . . . + z^, takes only three values on the set 
S'^ of vertices in S. In that case these three values are consecutive integers, 
say a + l,a,a — 1. We denote by S' (resp. S* resp. S°) the set of vertices 
where s takes the value a + 1 (resp. a resp. a — 1) and say that the vertices 
in S' are black, those in S° are white and those in S* are grey. Each of these 
sets is invariant under the translation action of L. 

6.10. Remark. I have at present no proof that for every rank 2 subgroup 
L C as in 12.11 the above construction indeed yields at least one perfect 
surface. On the other hand, in all examples I investigated the computer 
produced at least one perfect surface. It would also be very interesting to 
know which perfect surfaces can arise already in the first step 16.31 of the 
construction, by choosing A appropriately (instead of at random). 
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6.11. For every surface S in the list of 16.71 one has the graph F with set 
of vertices and arrows given by the diagonals in the squares, oriented 
such that in a square of type {i, j) the oriented diagonals are + and 
sign(det(bj, bj))(ei — ej). The advantage of having S embedded in is 
that the vertices and the arrows in these graphs are actual points and vectors 
in M^. If the surface S is perfect, every square in S has one vertex in S*, 
one in S° and two vertices in S*. One of its diagonals goes from the white to 
the black vertex and the other diagonal connects the two grey vertices. The 
graph F therefore is the disjoint union of two oriented graphs F'° and F*, 
with vertex sets S* U S° and S*, respectively. The graph F*° is a bi-partite 
graph, i.e. its vertex set is the disjoint union of two sets (the black resp. white 
vertices) and with edges connecting only vertices of different colors. 

There is a natural duality between F*° and F*: every vertex of F*° lies 
in a unique connected component of 5 \ F* and vice versa. Moreover the 
boundaries of these connected components are polygons and an arrow between 
two nodes in one graph is a common side of the corresponding polygons for 
the other graph. 

6.12. The perfect surface S is mapped homeomorphically onto Lj^ by the 
projection — ^ Lj^, which is dual to the inclusion L ^ Z^. This projection 
maps to bj, for i = 1, . . . , N. One may however pre-compose this projection 
with the linear map M^, nr-H^ii I < i < N. The composite 

II "^ i II 

linear map — >■ maps the tiling by squares on the surface S piecewise 
linearly and homeomorphically onto a tiling of by rhombi. It maps L 
isomorphically onto a lattice L in and it maps the graphs F*° and F* 
'isomorphically' onto L-periodic graphs F and F in the plane L^. In the 
literature F C is called a periodic dimer model or periodic brane tiling; 
see for instance [lEl [121 [131 [11 [15] . 

6.13. Definition. All structures in 16.111 are invariant under translation by 
vectors in L. Passing to the orbit space we obtain from a perfect surface S 
the 2-dimensional torus T = S/f^ (i.e. compact oriented surface of genus 1 
without boundary). Embedded in this torus are the two graphs F'° = F'yL 
and F* = ^J'L, which are dual to each other. Note that 16.121 also yields 
r = L^L, F«° = r°/L and F* = T/l. 

After Grothendieck one calls (T,F*,F*°) a dessin d'enfants fi8\ [22]. 

6.14. Example. Figure [TOl shows two drawings of the dessin (T, F*, F*°) for 
BiQ from Figure [3l or, rather, it shows a lifting of this dessin to the plane 
with basis so that the period lattice L becomes just Z^. It also shows F* as 
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an unembedded quiver. Note however that this is not the same as the quiver 
Q in case Biq in Figure\^ 

The quiver in Figure [10] is in fact the quiver for the singularity '^j'L^ 
where a generator of the cychc group acts on as multiphcation by the 
diagonal matrix diag(e^'^*/^, e^'^*/^, e^'^*/^). Also the brane tiling picture in 
Figure [TU] appears in [13] §3 in connection with the singularity C^^g. 

This difference between [13] and our approach comes, because we still have 
to perform the untwist, as is explained below. 




Figure 10: Tvjo versions of the dessin (T, F*, F*°) for Biq from Figure\^ The 
dashed square is the period parallelogram. The picture on the left shows how 
the quiver on the right is embedded in a 2-torus. 

6.15. The quiver F* in 16.131 is not the Pliicker quiver of L, defined in 12.51 
The latter appears in the current setting through the zigzag loops, as follows. 
Let 5 be a perfect surface. Fix i G {1, . . . , A^}. In every square on S which 
has two sides parallel to draw the line segment connecting the midpoints 
of these two sides. If the square has type this line segment is oriented 
so that it is a translate of the vector sign(det(bj, bj))ej. The union of these 
line segments projects to an oriented closed curve on T, called the i-th zigzag 
loop. Note that this zigzag loop may consist of several connected components. 

According to Proposition 16.51 the i-th and j-th zigzag loops on T intersect 
in exactly |det(bj,bj)| points. Thus we see that the nodes of the Pliicker 
quiver can be identified with the zigzag loops on T and the arrows between 
two nodes are the intersection points of the corresponding zigzag loops. The 
orientation of the arrows follows from the orientation of the zigzag loops and 
can in pictures be indicated as an over/undercrossing of the zigzag loops. 

Figure [11] shows the zigzag loops for the dessin in Figure [10] 

6.16. Final untwist. Pictures like Figure [TT] can be viewed as showing 
an oriented surface with boundary, in which the connected components of 
the zigzag loops are the connected components of the boundary and which 
is embedded in three space in a twisted way so that the black dots are on 
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Figure 11: The zigzag loops for the dessin in Figure fiOl 



one side of the surface and the white dots are on the other side. This surface 
comes with a tihng by hehces as shown in Figure [121 



b b b b 




W WW w 



Figure 12: Untwisting: from left to right: original square tile, helix tile, 
untwisted helix tile, untwisted helix completed to square tile. 

One can 'untwist' the hehx tiles and complete the untwisted helix tiles 
to squares as indicated in Figure [121 The effect is that the boundary cycles 
are being capped off by discs. The result is an oriented surface M. without 
boundary which is tiled by squares with one black, one white and two red 
vertices. The black-white diagonals give an embeding of the bipartite graph 
r*° into M.. The red-red diagonals give an embeding of the quiver Q into 
M. (if all zigzag loops have just one connected component; otherwise some 
points on the surface must be pinched together). 

The untwisting procedure is described in [3] §5 via examples, pictures and 
visual inspection. Here we want to build the 'untwisted surface' by simply 
reinterpreting the combinatorial data of the tiling of the torus T, i.e. the 
L-periodic tiling of the perfect surface S by unit squares taken modulo L (see 
16.91 16.13p . Let B (resp. W) denote the set of black (resp. white) vertices. 
There is a bijection between the set of tiles and the set E of arrows of the 
quiver Q and we can refer to tiles as e G i?. A tile e has one black vertex 
b(e) and one white vertex w(e). A tile e is the image of a unit square with 
sides parallel to two of the basis vectors ei, . . . , e^v; let us denote the indices 
of this (unordered) pair of basis vectors as a 2-element subset {r(e), r'(e)} of 
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{1, . . . , A^}. Thus the tihng on T yields the hst of quadruples 

M = {(b(e), w(e),r(e), r'(e))}ee£;. (26) 

Now, for every e G i? we take a unit square Dg and attach labels b(e), w(e), 
r(e), r'(e) and colors black, white, red, red, respectively, to its vertices such 
that the vertices labeled b(e) and w(e) are not adjacent. If two such squares 
□e and De' have equal labels on two adjacent vertices, we glue De and Dg' 
along the corresponding sides. The result of this glueing is a surface M. tiled 
with squares. ^ 

On the surface M. there is for every b G B one black point with label b 
and for every w G W there is one white point with label w. However for 
i G {1, . . . , A^} there may be several red points with label i. We finally obtain 
the desired surface M. by identifying for every i G {1, . . . , A^} the red points 
with label i. 

6.17. Proposition. The surface M., and hence also M., is oriented. 
Proof. Cut the square in four pieces like in the right-hand picture of 
Figure [121 Then for fixed b G B the small squares with black vertex b 
are glued together to a polygon isomorphic to the polygon in T formed by 
the unit squares with one vertex b. Similarly, for fixed w G W the small 
squares with white vertex w are glued together to a polygon isomorphic to 
the polygon in T formed by the unit squares with one vertex w but with the 
orientation reversed. For i G {1, . . . , A^} the small squares with red vertex i 
are glued together to a number of disjoint polygons, one for every connected 
component of the i-th zigzag loop and with oriented boundary 'equal to' that 
component. It is nowj^lear from Figure [T2] that these polygons are oriented 
consistently so that M. is oriented. ■ 

6.18. Conclusion. The surface Ai is tiled with squares so that the black- 
white diagonals form a bi-partite graph isomorphic to the bi-partite graph r*° 
in \6.13[ while the red-red diagonals form a graph isomorphic to the quiver Q. 



6.19. Example. The list M (126|) which our algorithm produces for model 
IV of dP^, i.e. example Biq from Figure [3l is 
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The surface with the corresponding tiling by kites (instead of squares) is 
shown in the left-hand picture in Figure [131 The right-hand picture in Figure 



30 



[T^ shows the bi-partite graph r*° and the quiver Q embedded in Ai. Note 
that this quiver is indeed the same as the one for example Biq in Figure [T31 




Figure 13: Surface M. with kite tiling (left) and quiver Q and bipartite graph 
V'° (right) for case Biq in Figure^ (= model IV of dP^). The surface M. is 
obtained by identifying opposite sides of the hexagon. 



7 Perfect matchings and the secondary fan. 

7.1. Definition, (see e.g. [T6l [20] ) A perfect matching P on a bi-partite 
graph r*° (also known as a dimer configuration) is a subset of the edges of 
r*° such that every node of r*° is incident to exactly one edge in P. 

7.2. Theorem. Consider a subgroup L C as in \2.1\ Let T*° be a bi- 
partite graph which is obtained from L C as explained in Section\^ and in 
Definition \6.13\ in particular. Let C be a 2-dimensional cone in the secondary 
fan of h G and let Lq be the corresponding collection of 2- element subsets 
of {1, . . . , N} as inlZEandlZM Identify E with the set of edges ofT'°. 
Then, with the notation as in ^2B) . the set 

Pc = {eeE\{r{e),r'{e)} G Lc } (27) 

is a perfect matching on r*°. 

Proof. Consider a perfect surface S C as in l6.9[ The map vr : — > Lj^ 
which is dual to the inclusion L ^ Z^ projects S homeomorphically onto the 
plane L^. The tiling of S by unit squares gives a tiling of Lj^ by parallelo- 
grams. A unit square with edges ej and e^ projects onto a parallelogram with 
edges bj and b^. Now let C be a 2-dimensional cone of the secondary fan 
(see 13. 2p and let w be a white vertex of the tiling. By translating the origin 
of the secondary fan to w one sees that there is exactly one parallelogram 
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with vertex w and with edges bj and hj such that C C M>obj + M>obj, 
i.e. such that {i,j} G Lq- This shows that there is exactly one element of 
Pc incident to w, namely the black- white diagonal of this parallelogram. A 
similar argument works for the black vertices. ■ 



7.3. Corollary. In the situation of Theorem 7.2 the number of black vertices 
and the number of white vertices are both equal to the number vol^ in [T^) . 
Proof. The number of black (resp. white) vertices of r*° is equal to the num- 
ber of edges in any perfect matching. This holds in particular for the perfect 
matching Pc corresponding to a 2-dimensional cone C in the secondary fan. 
The result now follows from fl27p . (jl]) and (IT^ . ■ 



The following result was also derived in [ij §5.1. 



7.4. Corollary. The genus of the surface M. in \6.1(^ is equal to the number 
of lattice points in the interior of the secondary polytope A in 
Proof. The number of 1-cells in the tiling oiA4 is twice the number of 2-cells. 
According to Proposition 16.51 the number of 2-cells is ^i<j<j<7v I det(bj, hj)\. 
The number of 0-cells is, according to Corollary 17. 31 and Formula (fT9!) equal to 
+ 2 J2{ij}eLc I bj)|. The desired result now follows from Corollary 

13.101 and a calculation of the Euler characteristic 2 — 2 genus ■ 



8 Kasteleyn matrix, bi-adjacency matrix, 
3-constellation and superpotential. 

8.1. In 16.121 we found the L-periodic bi-partite graph r*° in the plane L^. 
Periodic bi-partite graphs in the plane are also known as periodic dimer mod- 
els. One of the main tools for studying periodic dimer models is the Kaste- 
leyn matrix. It is usually defined as an adjacency matrix with extra factors 
X, x~^, y, y~^ for edges crossing the sides of a fundamental parallelogram and 
determined by visual inspection of a picture of the periodic dimer model; see 
e.g. [m [IS [15]. The purpose of the extra factors is to keep track of the ho- 
mology class of a closed path on the bi-partite graph r*° in the torus T (see 
I6.13p . We obtained P C as the projection of the black-white diagonals 
in the tiling of the perfect surface S by unit squares. The advantage of having 
S embedded in is that these diagonals are actual vectors in and that 
points lying on paths formed by these diagonals are simply given by vector 
addition. Thus all edges become "equally responsible" for the periodicity. An 
edge which is the diagonal in a square of type {i,j) contributes -|- e^. 

Thus we are led to a reformulation of the Kasteleyn matrix of the dimer 
model (T, P*°) in terms of the combinatorial data given as the list M in 



32 



Equation . As this same list is used in 16.161 for building the surface 
we call the same matrix also the bi-adjacency matrix of the dessin {A4, Q). 



8.2. Definition. From the list M = { (b(e), w(e), r(e), r'(e) ) }eeE in ([26]) 
and a function xv : E ^ C on the set of edges of the quiver Q we construct 
a matrix with rows indexed by the set of black vertices B, columns in- 
dexed by the set of white vertices W and with entries in the polynomial ring 
C[ui, . . . ,un\. 



(b,w)-entry of K"^ 



W{e) Ur(e)Ur'{e) ■ (28) 

e£E: b(e)=b, w(e)=w 



We use for the name Kasteleyn matrix of the dimer model (T, r*°) as well 
as the name bi-adjacency matrix of the dessin {Ai, Q). 



8.3. Example. Although the relation between the list M and the matrix 
is so straightforward, that there seems no point in writing here the matrix 
for Example I6.19[ we do present this here to compare it with results 
in the literature. 



bi-adjacency matrix 



for model IV of (iA 



ZU3U1U3 wjUiu^ 

w^uiu^ wiqu^uq W12U1UQ 

WijU^U^ VJQU1U4 ^gUi'U^ 

W^UiU^ Wi^U^Uq WuUiUq 

ZUgUiU^ tZ7i3tX2'U.5 rj72UiU2 

wiqUiu^ W1U1U2 'cui4U2Ue 

(with We denoting w(e)). In view of the discussion in Example 16.141 this 
bi-adjacency matrix for model IV of dPs should be somehow the same as the 
Kasteleyn matrix for the singularity C^^g in [13j Equation (3.5): 



K{z, w) 



Indeed, by dividing the first, third and fifth column of by and subse- 
quently setting M5 = , = w, ui = U2 = = uq = 1, and appropriately 
setting = =tl one obtains a matrix which up to a permutation of rows 
and columns is K{z, w). This way of getting K{z, w) from is nothing but 
passing from homogeneous to inhomogeneous coordinates. 
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8.4. Definition. (cf. [T8] Definition 1.1.1.) A 3 -constellation {E,ao,ai) 
consists of a finite set E and two permutations aQ,ai of E sucli tliat tlie 
permutation group generated by ctq, ai acts transitively on E. 

8.5. Notations. Tfie cycle notation {ii ^2 • ■ ■ ik-i ik) denotes the permuta- 
tion p such that p{ij) = ij^i for j = 1, . . . , — 1 and p{ik) = ii- By a cycle of 
a permutation a we mean an orbit of the group generated by a. We denote 
the set of cycles of a hj E^. In the product ar of two permutations a and r 
one first applies r ; so {aT){i) = cr(T(i)). 

8.6. Definition. For our purposes we do not need the most general notion 
of a superpotential for a quiver. The superpotentials we need are just no- 
tational reformulations of 3-constellations. The superpotential attached to a 
3- constellation {E,ao,(yi) is the following polynomial W in non-commuting 
variables X^. {e E E): 

where := Xj^Xj^ ■ . . . ■ Xj^. for a cyclic permutation 7 = (ii ^2 • • • "i/c)- 

8.7. Definition. The 3 -constellation {E,ao,ai) associated with the dessin 
{Ai, Q) in 16. 181 is the following. E is the set of arrows of the quiver Q. Every 
e G i? is a path on JH. Every connected component of \ Uei^EPe has an 
oriented boundary which can be viewed as a cyclic permutation of elements of 
E. Then ctq (resp. ai) is the composition of the cyclic permutations which are 
boundaries of connected components containing a black (resp. white) point. 

8.8. Remark. In the same way one associates a 3-constellation with the 
dessin (T, F*) in 16.131 It is obvious from the untwisting procedure in 16.161 
that, if {E,ao,ai) is the 3-constellation associated with the dessin {A4,Q), 
then {E, ao, CTi^) is the 3-constellation associated with the dessin (T, F*). 

8.9. It is obvious that the 3-constellation {E,aQ,ai) associated with the 
dessin {Ai,Q) contains the complete instructions for building Ai: for every 
cycle of (To and every cycle of o"i take a (convex planar) polygon with sides 
labeled by the elements in the cycle in their cyclic order. Next glue these 
polygons by identifying sides with the same label. The sides of these poly- 
gons correspond to arrows of the quiver Q and the vertices of these polygons 
correspond to vertices of Q. After the above procedure of glueing polygons 
along their sides one must still identify points which correspond to the same 
vertex of Q. The result is then J\4 with Q embedded in it. 
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8.10. The 3-constellation {E,aQ,(Ti) in 18.91 and the hst M in 16.161 both are 
completely determined by and do completely determine the dessin (A^,Q). 
One can, however, also describe the relation between {E,ao,ai) and M in a 
direct algebraic / combinatorial way. 

For the construction of {E,ao,ai) from M one first forms for b G B, 
w G W and z G {1, . . . , iV} the sets 

|o"b| = {e G -E I b(e) = b} , |crw| = {e E E \w{e) = w} , 
\zi\ = {ee E\ie {r{e),r'{e)}} . 

These sets have an unoriented cyclic structure: e, e' G |crb| (resp. e, e' G |o"w|) 
are neighbors if and only if {r(e),r'(e)} fl {r(e'), r'(e')} 7^ 0, while e, e' G |zj| 
are neighbors if and only if {b(e), w(e)} fl {b(e'), w(e')} 7^ 0. In order to put 
a consistent orientation on these cyclic sets we choose one of the two possible 
orientations of |zi|. For every b G B we have tt(|zi| fl \(Tb\) = or 2. In the 
latter case |zi|n|(Tb| consists of two elements, say e and e', which are neighbors 
both in |zi| and in |(Tb|. The orientation on |zi| then induces an orientation 
on IcTbl such that e is the successor of e' in |(Tb| if e is the successor of e' in |zi|. 
In this way the orientation on |zi| induces an orientation on every |(Tb| and 
every Icr^l which has a non-empty intersection with |zi|. Next this induces an 
orientation on every |zj| which has a non-empty intersection with any of the 
already oriented sets |o"b| or |crw|. And so on. It is because of the geometric 
background of M (i.e. the orientations shown in the right-hand picture in 
Figure [12]) that we can indeed go on. In the end all sets |(Tb| and \a^\ have 
an oriented cyclic structure and can be identified with cyclic permutations 
ab and cTw. The construction of the 3-constellation {E,ao,cri) is finished by 
setting 

beB wew 
the reason for inverting the white cycles is shown in Figure [T3] 



b 




w 



Figure 14: orientation of cycles in 3-constellation and orientation of Ai. 

To go from {E, ctq, cti) to M one views the cycles of the permutations ctq, 
(Ti and aoa^^ as elements of the sets E^^^, E„^ E^ ^-1, respectively, and then 
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sets M = {(b(e), w(e), r(e), r'(e))}egi? with 



b(e) = the cycle of 

w(e) = the cycle of 

r(e) = the cycle of 

r'(e) = r(ai(e)). 



ctq which contains e , 
o"i which contains e , 
aoai^ which contains e , 



8.11. Example. The procedure of 18.101 converts the list M in 16. 191 into the 
3-constellation {E,ao,ai) with 

E = {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18}, 
(To = (3,7,15)(4,12,16)(6,9,17)(5,11,18)(2,8,13)(1,10,14), 
ai = (4,8,15)(3,10,16)(5,7,17)(6,12,18)(1,9,13)(2,11,14). 

Since in this example the cycle zi = (1,10,3,7,5,11,2,8,4,12,6,9) meets 
each of the above 3-cycles the conversion went especially fast. 
The corresponding superpotential W is 

+ X4X12X1Q + X^XgXiT + X5X11X18 + X2X8X13 + X1X10X14 
— X4X8X15 — X3X10X16 — X5X7X17 — X6X12X18 — X1X9X13 — X2X11X14. 

9 The determinant of the bi-adjacency ma- 
trix and the secondary polygon 

9.1. Corollary 17.31 implies that the bi-adjacency matrix = (Kij) of the 
dessin {M., Q) is a square matrix. Its determinant is by definition 

detK"^ = y^^sign{r)Kir{i) 

T 

where r runs over the set of all bijections B ^ W and sign(r) is defined as the 
sign of the permutation Tq^t of B for some reference bijection tq : B ^ W. 
Of course, if we want to actually write IK™ as a matrix, we must choose 
bijections between B, W and the set of numbers {1, . . . , vol^} and that fixes 
tq. Changing the reference bijection multiplies detK™ by ±1, but that is for 
our purpose unimportant. Next note that the (b,w)-entry of K™ is non-zero 
if and only if there is an edge of r*° connecting the nodes b and w. So the 
only bijections r that contribute to the determinant of K™ are the perfect 
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matchings; see Definition 17. 1[ Thus we find 

detK"^ = ± ^ sign{P)Y[w{e)ur(e)Ur'{e) 

perfect matchings P egP 

= ± Yl sign(P)(n^(^))^^ (29) 

perfect matchings P \e€P / 

where 

P = ^(er(e)+er.(e)) G , (30) 

eeP 

and uP = ■ . . . ■ for p = (pi, . . . ,piy) E TL^ . 

9.2. Definition. The Newton polytope of a Laurent polynomial 



(fci,...,fciv)GZ^ 



is the polytope 

Newton(/) := convex hull ({(/ci, . . . , /cat) | C(fc,,...,fc^) 7^ }) 



9.3. Theorem. The Newton polytope of the determinant of the bi-adjacency 
matrix of the dessin (A^,Q) is the same as the secondary polygon of 
L C (see\Em): 

Newton(detK^) = S(L) . (31) 

The vertices of Newton(det K^) are the points Pc (see ^3^) ) given by the 
perfect matchings Pq with C a 2-dimensional cone in the secondary fan. 
Proof. If P and P' are two perfect matchings then the edges in P oriented 
black-to-white together with the edges in P' oriented white-to-black form a 
collection of closed loops on the torus T in 16. 131 The first homology group 
of this torus is L. Therefore P — P' lies in L C Z^. This shows that 
Newton (det K^) lies in a 2-dimensional plane in parallel to Ljr. 
From Equation (15U]1 one sees that for every perfect matching P 

N 



P = Y,^{e e P\i e {v{e),r'{e)}} 



1=1 



i.e. the z-th coordinate of the vector P equals the number of edges in the 
perfect matching P which intersect the i-th zigzag loop (cf. I6.15p . It fol- 
lows from Proposition 16.51 that the number of edges of the i-ih zigzag loop is 
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X^j=i I det(bj, bj)|. Since consecutive edges in a zigzag loop can never inter- 
sect the same perfect matching we see that for every perfect matching P and 
every i 

tt{e G P I z G {r(e),r'(e)}} < i E.li I ciet(b„b,)| . (32) 

Fix i and let C and C be the two 2-dimensional cones in the secondary fan 
containing the half-line M>obj. Then G Lc det(bj,bj) > and 

{i,j} G Lc -v^ det(bj,bj) < 0. Since X^jLi det(bj, bj) = we find 

J2 |det(b,,b,)| = Yl Idet(b„b,)| = 1 Ef=i I det(b„ b,)| . 

Thus we see that for i and the perfect matchings Pc and Pc Equation fl32l) 
is in fact an equality. 

Fix a 2-dimensional cone C in the secondary fan. It is bounded by two 
half-lines M>obj and M>obj/. The above argument now shows that Equation 
f l52]) is in fact an equality for the perfect matching Pc and i and for Pc and i' 
in place of i. Thus the point Pc is a vertex of Newton (det K^„). This together 
with (132|) proves: 

Newton(det IK^) = convex hull (^Pc \ C 2-dim. cone of secondary fan of . 

To complete the proof of fl3T|) we note that Equations ([5]), (|27|) and fl30l) 
together with Proposition 16.51 show Pn = ipc- ' 

9.4. Example. According to [S] §8 the quiver for model II of dP^ (i.e. case 
Bg in Figure Hj) affords two different superpotentials. This is confirmed by 
our algorithm, which yields the following two bi- adjacency matrices 

W2U1U2 ZUSU2U4 

^ ^uU^Uq UJ10U2UQ WgU2U^ 

^ WqUiUq + 'COJU2U2, Wi2,U4Uq W1U1U2 ZUuUsUi 

W2U1U2 W'^U2U^ W^UiU^ 

jj^j^ _ w^uiUi W1U1U2 ZUSU2U4 

^ UJ7U2U3 WuU^Uq W1QU2UQ Wi2U^U^ 

W11U3U4 WqUiUq ZUi^U^Uq W^lLiUs 

The corresponding superpotentials are: 

—-^^2-^7 -^12-^14-^^6 -'^8-'^13-'^10 — -^^1-^9-^5 — -'^3-'^ll-'^4 
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— — XiXgXi^X^ — XgXi^XiQ — X3X12X5 

After figuring out the rule for translating the edge labels one finds that the 
prepotentials in Eqs. (8.1) and (8.2) of [5] are Wa = Wi and Wb = W2. 

For both superpotentials the surface Ai has genus 1. Figures 18 and 19 
in [5] show pictures of the planar periodic bi-partite graph which is the lifting 
of r'° to the simply connected cover of A4. 

One can not expect the Kasteleyn matrices Ka and Kb in [S] Eq. (8.3) to 
be the same as the bi-adjacency matrices and K^, respectively, because 
of the "untwist" . and are Kasteleyn matrices of the dimer models 
before the untwisting (see 18.21) . Moreover Eqs. (8.4) and (8.5) of [5] show 
that the Newton polygons of the determinants of the Kasteleyn matrices Ka 
and Kb are different. On the contrary, for both and the determinant 
is a linear combination of the monomials 

^[3,3,1,0,0,1] ^[2,2,1,1,1,1] ^[1,1,1,2,2,1] ^[2,0,0,2,2,2] ^[3,1,0,1,1,2] ^[0,2,2,2,2,0] ^[1,3,2,1,1,0] 

So, in agreement with Theorem 19. 3[ the Newton polygons of det and 
det K2 coincide with the polygon in Figure H] case and 13. 7[ In spite of 
having the same Newton polygon det and det are not equal; indeed 
when all weights are 1 they differ in the coefficient of u'^'^'^'^'^'^l 



10 Bi-adjacency matrix with critical weights 
and the principal ^-determinant 

10.1. In [11] p. 297 Eq. (1.1) Gelfand, Kapranov and Zelevinsky define for a 
set A = (ai, . . . , a^v) C 'Z'''^^ as in 14.11 the principal A- determinant -E^(/^). 
The definition uses the Laurent polynomial 

TV 

U = J2u,^''\ (33) 

i=l 

where x™ = x^'^x^^ ■ . . . ■ x^J'^ for m = (mi, . . . , rrik+i) G Z'^'^^ and where 
the coefficients Ui, . . . ,un are variables. The name "principal ^-determinant" 
refers to the fact ([TTj p. 298 Prop. 1.1) that E^^ifA) can be written as the 
determinant of some exact complex, i.e. 

EaUa) = U^det M,)(~i)^ 
j 
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for certain matrices Mj. Another useful description of Ej\^{fj{) is given in [TT] 
p.299 Prop. 1.2: 

EMa)=± n ^AnrifAnrr^^^ (34) 

rcconvex hull(yi) 

where the product runs over all faces F of the primary polytope convex hull(^) 
(cf. 14.121) : m(r) is a multiplicity and A^nr(/Anr) is the {A HT) -discriminant 
of the Laurent polynomial /^nr = Xli aier The latter discriminant is a 

polynomial in the variables with aj G F. To define it (see [TT] p. 271) one 
needs the algebraic set V^nr which is the closure in C-^^'" of: 

{u G C-^"^ I 3x0 G (C*)'=+^ s.t. /^nr(xo) = ^^(xo) = 0, V*} . 

OXi 

Then, by definition, A_4nr(/^nr) = 1 if codimc-4nr(V^nr) > 1 and 

zero locus of A^nr(/^nr) = V^nr if codimc.Anr(V^nr) = 1 • 

So, Ejxi^fjx) gives the locus of the points (mi, . . . , ujy) G for which at least 
one of the Laurent polynomials /^nr has a critical point with critical value 0. 



10.2. Theorem. (fTT]p.302 Thm.1.4; cf.([I] 

The Newton polytope of Ej\^{fj) coincides with the secondary polytope S(^). 



10.3. Theorems 110.21 and 19.31 in combination with Formula fl20|) make one 
wonder whether for an appropriate choice of the weight zu the determinant of 
the bi-adjacency matrix is equal to the principal ^-determinant Eji^{fj{), 
up to the simple transformation necessitated by fl2U]) . In order to formulate 
this transformation we must make the dependence on lii, . . . ,ujv visible by 
writing W^{ui,..., un) and -E^(/^(mi, . . . , un))- 

After some experimenting with examples I found a very natural and simple 
choice for the weight that does the job. 

10.4. Definition. The critical weight for the arrows of the quiver Q is the 
function 

crit : E ^ Z>o , crit(e) = il{e' G E \ s{e') = s{e), t{e') = t{e)} (35) 



10.5. Conjecture. For every set A as in \4.1\ and every dessin {Ai, Q) (see 
\6.18\) constructed from A by the algorithm in Section the determinant of 
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the bi-adjacency matrix with critical weight and the principal A- determinant 
satisfy: 



[Ui- ...-UN 



)™'^ det W 



1 ) • • • ; 



EA{fA{Ul, - ■ ■ ,Un)) 



(36) 



10.6. Remark. In support of the above conjecture we can point out that the 
coefficients of the monomials corresponding to the vertices of the secondary 
polygons are, up to sign, the same for the two sides of Equation ^3^. Indeed, 
Theorem 1.4 of [lljp.302 gives the coefficient of the monomial in the principal 
^-determinant with corresponds to a vertex of the secondary polytope. Such 
a vertex corresponds to a maximal cone C in the secondary fan. In the 
notations of Definition 13.31 the formula in loc. cit. for the coefficient of the 
monomial corresponding to C reads: 



± n 



det(bi,bj^ 



dct(b,,bj)| 



(37) 



On the other hand. Theorems 17.21 and 19.31 show that the same vertex corre- 
sponds to a perfect matching Pq. From the role of perfect matchings in the 
computation of the determinant of the bi-adjacency matrix (see 19. ip one now 
easily checks that the coefficient of the monomial in the determinant, which 
corresponds to C, is (possibly up to sign) the same as dST 



In the remainder of this section we explicitly verify Conjecture 110.51 in 
some examples. 



10.7. Example. For case B2 



in Figure [3] the algorithm 



11-2 
-10 2-1 

in Section [6] yields in the following bi-adjacency matrix with critical weights: 



K' 



crit 



2UiU4 U1U3 
2U2U3 U2U4 
U1U2 + 3M3tt4 2ti2W3 2Wi«4 



One easily computes 

detK^"* = 27u[0'0'3'3l +4u[^'2'3'0] +4u[2'i'0'3] _ i8^[i,i,2,2] _ ^[2,2,1,1] 



Note that the exponents are the same as the coordinates of the vertices and 
the interior point of the secondary polygon in 13.51 For the computation of 
the principal ^-determinant we note that in this case (see also Example 14.61) 



A = (ai, a2, as, a4) = 





1 




1 




1 




1 


( 





1 


3 


1 


1 


1 


2 
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and hence = Xi{ui + U2xl + u^X2 + u^xl). The primary polytope is shown 
in Figure [H The two boundary points of this primary polytope contribute 
factors Ui and U2 to the principal ^-determinant (cf. Equation flM|) ). The 
contribution from the full primary polytope is the very classical discriminant 
of the cubic polynomial ui + U3X2 + U4XI + U2xl, which is (e.g. |TT] p. 405) 

21UiU2 + AUiU^ + Au^U2 — "^3^4 — \'SUiU^U4U2 ■ 

Thus we find 

EAifA) = 27u[=^'3'0'0l + 4uP'i.M + 4u[i.2,3,o] _ ^[1,1,2,2] _ l^^%2,l^] 
This shows that conjecture \10.5\ holds in this case. 

10.8. Example. When critical weights are used the two bi-adjacency ma- 
trices in Example 19.41 become 

2uiU2 U2U4 U1U4 

U5UQ U2UQ U2U^ 

UiUq + U2U3 2UiU2 U3U4 

UsU5 

2uiU2 U2U5 
U1U4 2uiU2 U2U4 
U2U3 u^ue U2U6 U3U5 

UiUq U4UQ U1U3 

One easily computes det K^"^'* and det Kf^^ and finds that both are equal to 

^^[3,3,1,0,0,1] _^ 2u[^'^'^'^'^'^l + u[^'^'^'^'^'°l + u[^'^'°'^'^'2l - 
_^[2,0,0,2,2,2] _ ^[0,2,2,2,2,0] _ g^[2,2,l,l,l,l] 

This gives a remarkable contrast with the last line of Example 19.41 and un- 
derlines the role of the critical weight in making Conjecture 110.51 reasonable 
for every dessin obtained from A. 

Transforming the above determinant as in the left-hand side of Equation 
(IMD yields 

^^[1,1,3,4,4,3] _^ 2u[^'^'^'^'^'^l + u[3'^'2'3'3'4] _^ ^[1,3,4,3,3,2] _ ^[2,4,4,2,2,2] _ 

_^[4,2,2,2,2,4] _ g^[2,2,3,3,3,3] 

= UiU2ulululul{u4U5 — UiU2){4:U3U4U5Uq + ulul — 2uiU2UsUq + ulul). (38) 

From the matrix Bg in Figure H] one easily finds A and the corresponding 
polynomial 

Ia = "Uia^i + U2X1X2 + U3X1X3 + u^XiXi + U5XiX2X'^^ + UGX1X2X3 . 
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The primary polytope convex hull(^) is 3- dimensional and has seven 2-di- 
mensional faces, six of which are triangles and one is a quadrangle. Its only 
integer points are its vertices and these correspond to the monomials of fj^. 
One can compute the discriminants for the polynomials supported by the faces 
of the primary polytope and multiply these to get the principal ^-determinant 
EaUa) ill Equation flMl) . The result is exactly as in fl55]) . with the 4-term 
factor coming from the primary polytope itself, the 2-term factor coming from 
the quadrangle face and the other factors with multiplicities coming from the 
vertices. This shows that conjecture \10.5\ holds in this case. 



10.9. Further examples. Exactly as in Example I1U.7I one can verify Con- 
jecture dOSl for other "four-nomials" of degree < 5 (see Example 14. 7p by using 
the formulas in [llj for discriminants of quartic and quintic polynomials. 

The method of Example 110.81 can be successfully employed to verify Con- 
jecture [iniS] for, for instance, cases B^, Bq, Biq in Figure HI It also works 



for B 



and B 



-311100 
-300011 
closely related to Biq and Bg, respectively. 




1 



1 




-2 
1 





-2 



which are 
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